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Abstract 

Exactness of the lowest order supersymmetric WKB (SWKB) quantization 
condition f^'^ \/E — u>'^{x)dx = nhir, for certain potentials, is examined, using 
complex integration technique. Comparison of the above scheme with a sim- 
ilar, but exact quantization condition, §^p[x,E)dx = 2'Knfi, originating from 
the quantum Hamilton- Jacobi formalism reveals that, the locations and the 
residues of the poles that contribute to these integrals match identically, for 
both of these cases. As these poles completely determine the eigenvalues in 
these two cases, the exactness of the SWKB for these potentials is accounted 
for. Three non-exact cases are also analysed; the origin of this non-exactness 
is shown to be due the presence of additional singularities in \J E — a;^(x), 
like branch cuts in the x— plane. 
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Exactness of the lowest order supersymmetric WKB (SWKB) approximation scheme for 
a class of potentials has been the subject of wide discussions in the literature [1]. Explicit 
calculations, akin to that of the celebrated harmonic oscillator case [2], first revealed the 
absence of 0(^^) to 0(^^) corrections to the lowest order result [3]. Subsequently, for one 
of the above potentials, it was shown that all higher order corrections vanish identically [4] . 
A discrete reparameterisation symmetry, the so-called shape invariance, of these potentials, 
has been thought to be responsible for this exactness [5]. Recently, a new class of shape 
invariant potentials have been found, for which the exactness of SWKB quantization does 
not hold true [6]. In this light and considering the importance of this quantization scheme 
in the potential problems, a deeper analysis of the origin of this exactness is of utmost sig- 
nificance. 

In supersymmetric (SUSY) quantum mechanics [7] , a pair of Hamiltonians, H± with po- 
tentials V±{x), given by uj'^{x) ± hduj{x)/dx, (2m = 1), have identical bound state spectra, 
except for the ground state; u!{x) here, is the superpotential in SUSY quantum mechanics. 
In the case of unbroken SUSY, conventionally, the energy of the ground state of which 
is unpaired, is taken to be zero. It was observed in [1], that a semi-classical WKB type 



approximation J^^ ^ E — (jj'^(x)dx = nhir, making use of {V-{x) + V-^-{x))/2 = u!'^{x), gave 
exact energy eigenvalues, for many well known potentials. This is to be contrasted with 
the conventional WKB scheme, which reproduced the exact spectra only for the harmonic 
oscillator and the Morse potentials. 

In what follows, we will first use a complex integration technique, to evaluate the SWKB 
integral given by. 



Here, xi and X2 {xi < X2) are the turning points, which are the real roots oi E — uj'^{x) = 
and C is a counter clockwise contour, enclosing the branch cut between xi and X2- Note 
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JswKB = ~ y E ~ 'jj'^{x)dx = — f y E — uj'^{x)dx = nh 




(1) 



that the integrand is a double valued function, which is taken to be positive just below the 
cut, joining the two turning points [8]. The above integral will, then, be contrasted with an 
exact quantization condition. 



originating from the quantum Hamilton- Jacobi (H-J) formalism [9,10]. Here, J{E) is the 
quantum action variable and p{x, E) is the quantum momentum function (QMF) and C 
is a counter clockwise contour in the complex x— plane, enclosing the real line between the 
classical tuning points. The classical turning points are the real values of x, for which 
E — V{x) = pI{x,E) vanishes, and encloses the region between which the classical motion 
takes place. In a previous paper [11] we used the quantization condition given by Eq. 
to obtain the energy eigenvalues for a set of SUSY potentials. This was done by rewriting 
the integral J{E) in Eq. as sum of integrals over contours enclosing the singular points 
of QMF. 

In this paper we evaluate Jswkb by complex integration technique in a manner similar 
to that used in [11] for J{E). Interestingly, we find that the locations and the residues of the 
poles of QMF, which determine the eigenspectra for the exact quantization scheme match 
identically with the poles appearing in the evaluation of the integral in Eq. (|^. This, then, 
accounts for the SWKB scheme being exact in these cases. 

In this paper we also analyse three potentials for which the energy eigenvalues are known, 
but SWKB fails to give the exact answer [5,12]; It is seen that the presence of additional 
singularities, like branch cuts, prevents an exact calculation of SWKB integrals for all these 
cases. Intriguingly, in the above examples, it is found that, if one ignores the contributions of 
some of the branch cuts in the non-classical region, the exact answer is reproduced. Hence, 
for these potentials, the difference between the exact answer and the SWKB approximation 
scheme, is quantified; this may be quite useful in developing new approximation schemes for 




(2) 
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non-exactly solvable problems. 

The QMF satisfies the quantum H-J equation (2m = 1) 



I OX 

^pI{x,E) . (3) 
Here, V{x) is any given potential and Pc{x, E) is the classical momentum function, defined 



so as to have the value +y E — V{x) just below the branch cut which joins the two turning 
points. QMF p{x, E) obeys 

limp(a;, i?) — > Pc(a;, -E") . (4) 
?i— >o 

This can be thought of as a manifestation of the correspondence principle or as a boundary 
condition on the QMF. The quantum H-J equation is equivalent to the Schrodinger equation, 
under the substitution, 

^1 dil}{x,E) 
^^"'^^ = 7 ^(a;,E) dx ■ 

The nodes of the wave function E) are known to appear between the classical turning 

points. QMF E) will then have poles at the corresponding points, hence J{E) which 

counts this number, will obey the exact quantization condition, given in Eq. (Q). If we 

substitute an % expansion for E), the quantization condition (j^) coincides with the well 

known Dunham's formula. 

We explicitly workout the SWKB integral in Eq. (|I|) for a set of trigonometric and hy- 
perbolic potentials and tabulate the results in Table 1. In what follows, the steps required to 
solve one of these potentials are outlined. We then compare these results with our previous 
calculation [11] using the quantum action variable. 

We shall now evaluate Jswkb for the Eckart potential, 

VI (x) = + — + A{A - ah)cosech^ax - 2B coth ax . (6) 



whose superpotential is, 

uj{x) = -Acothax + B/A , B > A^ and < z < oo . (7) 

We use the mapping y = exp(ax) in which case, 

V( ^ A^M^\ 4A(A - ah)y' 2B{y^ + 1) 



and 
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The quantization condition in Eq. (|l]) for the SWKB approximation becomes, 



1 r JE-iu^{y) 

JswKB = 7, — f dy = nh . (10) 

zvra Jci y 

It is clear from Eq. (|T0]) that, the integrand has singularities at y = and y = ±1. The 
roots of E — uj'^{y) = 0, give the branch points; they are four in number, two of which lie in 
the classical and the other two in the non-classical {y < 0) region of the complex ?/— plane. 
To evaluate the above integral, consider a counter clockwise contour integral Jr^ for a circle 
of radius R, which is such that it encloses all the singular points of the above integrand 
(see Fig. 1). Then, 

= JsWKB + Jc2 + J-n + -^72 + ^73 • (11) 

Here, is the integral along the counter clockwise contour C2 enclosing the branch cut in 
the Rey < region. J^^, and are the integrals along contours 71, 72 and 73 enclosing 
the singular points at y = 1, ?/ = — 1 and y = respectively. It may be noticed that the 
symmetry y —y of E — uj'^{y) implies, 

JswKB = Jc2 ■ (12) 

For the pole at ?/ = 0, one gets. 



iJE ~ A^ - B'^/A^ -2B 
J,, = ^ ^ . (13) 



Similarly for poles at y = ±1, 



J-y2 = J-y3 = — . (14) 



Now for the calculation of Jr^ one more change of variable z — 1/y is sought; so the 
singularity at y ^ oo is mapped to the singularity at z = 0. The contour integral for the 
pole at ^ = is, 



-i^E-A^-B'^/A^ + 2B 
Jtr = ■ (15) 



a 



The SWKB quantization condition 



^r«-EJ> = 2n^ , (16) 
p=i 



when inverted for E gives 



E. = A'^BVA^-^-^^-(ru,n + Af . (17) 

The calculation of eigenvalues for other SUSY potentials proceeds along similar lines and 
the results are summarised in Table I. 

Comparing the above computation of Jswkb, with that of J{E) in our previous paper, we 
find that the singularity structure of \JE — uS^ other than the branch cuts, match exactly 
with that of the fixed poles of p(x, E) in the quantum H-J formalism [11]. In all the cases 
tabulated in Table I, the location of the fixed poles and the values of the corresponding 
residues are identical for both the SWKB approximation and the quantum H-J formalism. 
The contour integral J(-E'), matches exactly with the Jswkb- This then explains the rather 
mysterious result of the exactness of SWKB approximation for the potentials in Table I. 



We next consider the cases where lowest order SWKB fails to give exact answer. We at 
first look at the potential [12] 
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with 

, , 2x^ + 3x^ -x^ -6 

The SWKB quantization condition reads 

JswKB ^ ^ E - uj^{x)dx = nh , (20) 

It is noticed that the above integrand becomes zero at twelve points in the complex plane 
giving rise to six branch cuts. One of these branch cuts is included in the contour C. The 
integral around a branch cut on the negative axis is equal to the above integral due to the 
symmetry x — > —x. 

The integrals around the other four branch cuts cannot be related to the integral Jc- We 
use JoBC to denote the sum of the integrals along contours enclosing the four remaining 
branch cuts. 



We now try to proceed, as far as possible, parallel to the Eckart potential. Let Jr^ be a 
circular contour with the center at origin and radius large enough to enclose all the singular 
points of the integrand. Note that the integrand has poles a,t y — 0; y — ±i and y — ±V2i. 
The contour integral Jr^ is given by 

Jpj} = 2nh + Jqbc + <^7i + ~^ ~^ J-1A ~^ ' (2-'-) 

or 
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■^Fh - E >^7p = 2n/i - JoBC , (22) 

are the corresponding contour integrals for the contours 
111 72, 73, 74 and 75, around y = 0; +i; — i; and - -\/2i 



Various contour integrals can be evaluated as before and we get J^^ = 3/2 , J^^ = ~1 , "As ~ 
—1,J^^ — 1 and J-yg = 1. The signs of the residues needed have been fixed as for the Eckart 



potential. Jr^ is again calculated making the mapping y ^ 1/z and looking for the residue 
of the pole at 2; 0. The corresponding contour integral is found to be, 

Jr^ = 2E + 3/2 . (23) 

The SWKB quantization condition is equivalent to Eq. ( |2^). and this condition involves 
the unknown integral Jqbc- Therefore, this complex variables approach cannot be used to 
complete the calculation of energy eigenvalues in SWKB approximation for the potential 
( |18|) . However, it is easily checked that if we substitute the values of Jr^, J^^, J-yai Jn 
and in 

Jv^-jZJ,. = '^nn , (24) 
p=i 

we obtain the exact eigenvalues. 

Therefore, it is clear that if the sum of the contributions coming from other branch cuts is 
nonzero, the SWKB approximation will not be exact. 



A similar analysis can be repeated for some of the other potentials. It is known that [12,5] 
for 



and 

with super potentials. 



-XMy + 1) + j(l - A')|2 - (7 - \')y' + 5(1 - A=)/] 



uj{x) 



and 



- IQx^ - 56x3 _ ^Q8^2 _ 240a; - 192 
4x(a;2 + 2x + 2){x^ + Gx^ + 16a; + 24) 



uj{x) = -{l-\^)y{y^-l)^i^\^y 



(25) 



(26) 



(27) 



(28) 



SWKB does not give exact answers. We once again see that SWKB is not exact because of 
the contributions from the other branch cuts, that are present in both these cases. 
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V. CONCLUSION: 

To conclude, we have shown using the complex contour integration technique, for several 
solvable potentials in one dimension, that SWKB gave exact results. We then compared it 
with the answers gotten using the quantum H-J method, from the exact quantization condi- 
tion. This has provided some useful insight: It is the similarity of the singularity structure 
of p{x, E) and \JE — lJ^ in the non-classical regions of the a;— plane, namely the matching 
of the poles and the residues that is responsible for this exactness. Table I contains all the 
relevant details about the steps involved in the calculation. 

For the non-exact potentials, it is the inability to deform the contour appropriately because 
of the presence of other poles and branch cuts in the complex x— plane that prevents an 
exact solution in SWKB approach. 

Acknowledgements: We acknowledge useful discussions with Prof. V. Srinivasan and Dr. 
C.N. Kumar. 
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Table I : Hyperbolic and trigonometric potentials. The mapping used for hyperbolic 
potentials is y = exp(Q;a;), while for the trigonometric ones y = exp{iax) is being used; 
a is real and positive. 
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Figure Captions 

1. Fig.l. Contour for the Eckart potential problem. 



